FINI
/CLE
/FILNAME,e422
/TITLE,EATHQUAKE ANALYSIS OF SPHERICAL TANK

!βλMM
S_R=15700/2          		!ǰ뾶
Z_R=560/2            		!֧ܰ뾶
D_R=986/2            		!֧װ뾶
ZCB_R=3/2*Z_R        		!Ӱ뾶
UPZ_H=2700           		!֧
LOWZ_H=10340         		!֧λ
Q1_H=2/5*UPZ_H       		!ǿIֱλ

!ϲ
S_E=2.094E5          		!ģ
S_U=0.262            		!ǲɱ
S_DENS=7.83E-6       		!ܶ

Z_E=2E5              		!ģ
Z_U=0.3              		!ϲɱ
Z_DENS=7.8E-6        		!ܶ

!ʵ
S_T=44               		!ǡӰ
D_T=54               		!֧װ
Z_T=10               		!
TKWALL=10            		!֧ܺ
LG_R=50/2            		!˰뾶

/PREP7
ET,1,93              		!SHELL93
ET,2,16              		!PIPE16
ET,3,8               		!LINK8

KEYOPT,2,6,2

MP,EX,1,S_E          		!ǲģ
MP,PRXY,1,S_U        		!ǲϲɱ
MP,DENS,1,S_DENS     		!ǲܶ

MP,EX,2,Z_E          		!ģ
MP,PRXY,2,Z_U        		!ϲɱ
MP,DENS,2,Z_DENS     		!ܶ

R,1,S_T              		!ǡӰ
R,2,D_T              		!֧װ
R,3,Z_T              		!
R,4,Z_R,TKWALL       		!֧
R,5,3.14*LG_R**2     	        !

!-------------------
!----νģ-------
!-------------------
!-------ǽģ
K,1,S_R              		!ؼ1(S_R,0,0)
K,2,,-S_R            		!ؼ2(0,-S_R,0)
K,3,,S_R            		!ؼ3(0,S_R,0)
CSYS,1               		!
L,2,1		     		!ӹؼ2,1Բ1
L,1,3		     		!ӹؼ1,3Բ2
AROTAT,1,2,,,,,2,3,18		!1,22,3Źؼᣬת18Բ

!-------ǰ彨ģ
WPOFFS,S_R			!Xƫƹƽ棬S_R
CSWPLA,11,1			!ڹƽԭ㴦岢11žֲϵ
K,5,Z_R				!ؼ5
K,6,Z_R,90			!ؼ6
L,5,6				!ӹؼ5,6Բ
AROTAT,6,,,,,,1,6,180,1		!61,6Źؼᣬת180Բ

APTN,2,3			!23滥

!-------а彨ģ
K,100,S_R/3,180,-Z_R		!ؼ100
K,101,,,-Z_R			!ؼ101
L,100,101			!ӹؼ100,101

WPROTA,,90			!Xתƽ棬90
CSWPLA,12,1			!ڹƽԭ㴦岢12žֲϵ

K,102,Z_R			!ؼ102
L,102,101			!ӹؼ102,101
LCOMB,2,7			!ϲ27

K,1000,,,UPZ_H			!ؼ1000
K,1001,,,LOWZ_H			!ؼ1001
L,1,1000			!ӹؼ1,1000
L,1000,1001			!ӹؼ1000,1001
ADRAG,2,,,,,,7			!7߶·2ߣ

WPOFFS,,,Q1_H			!Zƫƹƽ棬Q1_H
ASBW,ALL			!ùƽָǰѡ
WPOFFS,,,UPZ_H-Q1_H		!Zƫƹƽ棬Q1_HUPZ_H-Q1_H
ASBW,ALL			!ùƽָǰѡ

APTN,8,9			!89滥
APTN,2,10			!210滥

!-----֧װ
CSYS,12				!12žֲϵ
K,1002,D_R,,LOWZ_H		!ؼ1002
L,1001,1002			!ӹؼ1001,1002
AROTAT,18,,,,,,1,6,180,2	!181,6Źؼᣬת180Բ

!-----ɾ
CSYS,2				!2ϵ
ASEL,S,LOC,X,0,S_R-1		!ѡǰ뾶S_Rڵ
ADEL,ALL,,,1			!ɾж

CSYS				!ֱϵ
ALLSEL				!ѡԪ

NUMMRG,ALL			!ϲԪ
NUMCMP,ALL			!ѹԪر

!----ǿI
A,12,15,10			!ʹ121510Źؼ㣬
!----Ӱ
A,9,17,13,20			!ʹ9171320Źؼ㣬

!-----ָ棬㻮
WPCSYS,,0			!ƽѿϵ
WPROTA,,,4			!Yתƽ棬4
ASBW,ALL			!ùƽָǰѡ

WPCSYS,,0			!ƽѿϵ
WPROTA,,-90			!Xתƽ棬-90
WPROTA,,,-4			!Yתƽ棬-4
ASBW,ALL			!ùƽָǰѡ

WPROTA,,,28			!Yתƽ棬28
ASBW,ALL			!ùƽָǰѡ

WPROTA,,,66-10			!Yתƽ棬56
ASBW,ALL			!ùƽָǰѡ

WPROTA,,,2*10-180		!Yתƽ棬-160
ASBW,ALL			!ùƽָǰѡ

ALLSEL				!ѡԪ
AGLUE,ALL			!ճ

NUMMRG,ALL 			!ϲԪ
NUMCMP,ALL			!ѹԪر

!-------------------
!-----------
!-------------------
MSHAPE,0			!ָԪ״Ϊı
MSHKEY,1			!ָΪӳʽ
!Ӱ
ASEL,S,AREA,,10			!ѡ10
AESIZE,ALL,200			!趨ѡ浥ԪߴΪ200
TYPE,1				!ʹ1ŵԪ
MAT,2				!ʹ2Ų
REAL,1				!ʹ1ʵ
LCCAT,35,36			!35,36
AMESH,ALL			!ѡ
ALLSEL				!ѡԪ
!֧װ
ASEL,S,AREA,,2,6,4		!ѡ2,6
AESIZE,ALL,200			!趨ѡ浥ԪߴΪ200
TYPE,1				!ʹ1ŵԪ
MAT,2				!ʹ2Ų
REAL,2				!ʹ2ʵ
AMESH,ALL			!ѡ
ALLSEL				!ѡԪ
!
CSYS,2				!2ϵ
ASEL,S,LOC,X,S_R+1,2*S_R	!ѡ뾶ΪS_R+12*S_R֮
ASEL,U,AREA,,10			!ӵǰѡȥ10
ASEL,U,AREA,,2,6,4		!ӵǰѡȥ2,6
AESIZE,ALL,200			!趨ѡ浥ԪߴΪ200
TYPE,1				!ʹ1ŵԪ
MAT,2				!ʹ2Ų
REAL,3				!ʹ3ʵ
LCCAT,64,63			!64,63
LCCAT,14,6			!14,6
LCCAT,63,65			!63,65
AMESH,ALL			!ѡ
ALLSEL				!ѡԪ
!
ASEL,S,LOC,X,S_R		!ѡ뾶ΪS_Rϵ
ASEL,R,LOC,Z,-4,0		!ӵǰѡѡFIǶ-40֮
ASEL,R,LOC,Y,-30,4		!ӵǰѡѡTHETAǶ-304֮
AESIZE,ALL,200			!趨ѡ浥ԪߴΪ200
TYPE,1				!ʹ1ŵԪ
MAT,1				!ʹ1Ų
REAL,1				!ʹ1ʵ
MSHAPE,0			!ָԪ״Ϊı
MSHKEY,0			!ָΪʽ
AMESH,ALL			!ѡ
ALLSEL				!ѡԪ
!
ASEL,S,LOC,X,S_R		!ѡ뾶ΪS_Rϵ
ASEL,R,LOC,Z,-18,-4		!ӵǰѡѡFIǶ-18-4֮
ASEL,R,LOC,Y,-30,4		!ӵǰѡѡTHETAǶ-304֮
MSHAPE,0			!ָԪ״Ϊı
MSHKEY,1			!ָΪӳʽ
!γ߷򻮷
LESIZE,45,,,3			!ָ45߻Ϊ3
LESIZE,22,,,6			!ָ22߻Ϊ6
LESIZE,24,,,9			!ָ24߻Ϊ9
LESIZE,51,,,3			!ָ51߻Ϊ3
LSEL,S,LINE,,27			!ѡȡ27
LSEL,A,LINE,,39,41,1		!ǰѡм39,40,41
LSEL,A,LINE,,49			!ǰѡм49
LESIZE,ALL,,,2           	!ָѡ߻Ϊ2
AMESH,ALL			!ѡ
ALLSEL				!ѡԪ
!´
ASEL,S,LOC,Y,4,80		!ӵǰѡѡTHETAǶ480֮
LSEL,S,LINE,,59,61,1		!ѡȡ59,60,61
LESIZE,ALL,,,10,0.3		!ָѡ߻Ϊ10,Ϊ0.3
AMESH,ALL			!ѡ
ALLSEL				!ѡԪ
!´
ASEL,S,LOC,Y,-80,-30		!ӵǰѡѡTHETAǶ-80-30֮
ASEL,R,LOC,X,S_R		!ӵǰѡѡ뾶ΪS_Rϵ
LSEL,S,LINE,,54,56,1		!ѡȡ54,55,56
LESIZE,ALL,,,10,0.5		!ָѡ߻Ϊ10,Ϊ0.5
AMESH,ALL			!ѡ
!
ASEL,S,LOC,Y,80,90		!ӵǰѡѡTHETAǶ8090֮
ASEL,A,LOC,Y,-90,-80		!ǰѡмTHETAǶ-90-80֮
MSHAPE,0			!ָԪ״Ϊı
MSHKEY,0			!ָΪʽ
AESIZE,ALL,200			!趨ѡ浥ԪߴΪ200
AMESH,ALL			!ѡ
ALLSEL				!ѡԪ

!ɾ
LSEL,S,LCCA			!ѡ
LDEL,ALL			!ɾǰѡ
ALLSEL				!ѡԪ

!X-Yƽ澵
CSYS				!ֱϵ
ARSYM,Z,ALL			!X-Yƽ澵Ƶǰ棬ͬʱ

!֧
LSEL,S,LINE,,8			!ѡȡ8
LESIZE,ALL,,,10			!ָѡ߻Ϊ10
TYPE,2				!ʹ2ŵԪ
MAT,2				!ʹ2Ų
REAL,4				!ʹ4ʵ
LMESH,ALL			!ѡ
ALLSEL				!ѡԪ

NUMMRG,ALL 			!ϲԪ
NUMCMP,ALL			!ѹԪر

!ʾ1/10ģ
/VIEW,1,1,1,1			!ָһ
EPLOT				!ʾԪ

!10γȫģ
CSYS,5				!YΪתϵ
AGEN,10,ALL,,,,36		!THETA36ΪƫƵǰ10()ͬʱ
LGEN,10,8,,,,36			!THETA36Ϊƫ810()ͬʱ
ALLSEL				!ѡԪ
NUMMRG,ALL 			!ϲԪ
NUMCMP,ALL			!ѹԪر

!˵Ԫ
TYPE,3				!ʹ3ŵԪ
MAT,2				!ʹ2Ų
REAL,5				!ʹ5ʵ
E,3185,30710			!3185,30710Žڵ㣬ɵԪ
E,3193,30702			!3193,30702Žڵ㣬ɵԪ
E,30702,30719			!30702,30719Žڵ㣬ɵԪ
E,30710,30711			!30710,30711Žڵ㣬ɵԪ
EGEN,8,9,10583,10584, 1,,,,,,,36!THETA36ΪƫӦڵΪ910583,10584ŵԪ8()
E,3185,30782			!3185,30782Žڵ㣬ɵԪ
E,3193,30774			!3193,30774Žڵ㣬ɵԪ
CSYS				!ֱϵ

!ɲģ
ALLSEL				!ѡԪ
NUMMRG,ALL 			!ϲԪ
NUMCMP,ALL			!ѹԪر
/VIEW,1,1,1,1			!ָһ
EPLOT				!ʾԪ
SAVE				!

!-------------------
!-----------
!-------------------
!֧ʩӹ̶Լ
ASEL,S,LOC,Y,-LOWZ_H		!ѡȡYΪ-LOWZ_H
DA,ALL,ALL			!Լɶ

!غ
*DIM,TJX,ARRAY,2,190,0		!2X190TJX洢ˮƽ
*DIM,TJY,ARRAY,2,190,0		!2X190TJY洢ֱ

*CREATE,ANSUITMP		!ANSYSʱļANSUITMP
*VREAD,TJX(1,1),'ACELX','TXT',' ',190	!ACELX.TXTݵTJX
(E9.3,E11.3)				!ݸʽ
*END				!ļANSUITMP
/INPUT,ANSUITMP 		!ӺļANSUITMPж

*CREATE,ANSUITMP		!ANSYSʱļANSUITMP
*VREAD,TJY(1,1),'ACELY','TXT',' ',190	!!ACELY.TXTݵTJY
(E9.3,E11.3)				!ݸʽ
*END				!ļANSUITMP
/INPUT,ANSUITMP 		!ӺļANSUITMPж

/SOLU		
ANTYPE,4			!ָΪ˲̬ѧ
TRNOPT,FULL			!˲̬ѧFULL
OUTRES,BASIC,ALL		!ÿһ

!ʩӵغ

*DO,T,1,190,1			!ʼt1190ѭ
TIME,0.1*T			!趨غɲĽʱΪ0.1*T
KBC,0				!ָغΪغ
NSUB,1				!ָغӲΪ1
ALPHAD,0.05 			!趨ALPHAΪ0.05
BETAD,0.01			!趨նBETAΪ0.01

ACEL,TJX(2,T),TJY(2,T)
ALLSEL				!ѡԪ
SOLVE				!һغɲ
*ENDDO				!ѭ
SAVE				!

!-------------------
!----ͨú-----
!-------------------
FINI
/POST1
SET,,,,,12			!ӽļ12ʱ̵
PLNSOL, U,SUM, 0,1.0     	!ʾλͼ
PLNSOL, U,X, 0,1.0	     	!ʾXλͼ
PLNSOL, U,Y, 0,1.0	     	!ʾYλͼ

!ʾʱȥӦеӰ嵥Ԫ
ESEL,S,MAT,,2			!ѡϺΪ2ĵԪ
ESEL,R,REAL,,1			!ӵǰѡѡʵΪ1ĵԪ
ESEL,INVE			!ѡǰԪ			
PLNSOL, S,EQV, 0,1.0    	!ʾEQUIVALENT STRESS

!ǲֺ
ESEL,S,MAT,,1            	!ѡǵԪ
PLNSOL, S,EQV, 0,1.0     	!ʾEQUIVALENT STRESS
ETABLE,Q_MX, SMISC, 4    	!XصԪ
ETABLE,Q_MY, SMISC, 5    	!Y
ETABLE,Q_TX, SMISC, 1    	!X
ETABLE,Q_TY, SMISC, 2    	!Y

PLETAB,Q_MX,NOAV         	!ʾXͼ
PLETAB,Q_MY,NOAV         	!ʾYͼ
PLETAB,Q_TX              	!ʾXͼ
PLETAB,Q_TY              	!ʾYͼ

!а岿ֺ(Ӱ)
ESEL,S,ENAME,,93		!ѡԪΪSHELL93ĵԪ
ESEL,U,MAT,,1			!ӵǰѡȥϺΪ1ĵԪ
ESEL,U,REAL,,2			!ӵǰѡȥʵΪ2ĵԪ
ESEL,U,REAL,,1			!ӵǰѡȥʵΪ1ĵԪ

PLNSOL, U,SUM, 0,1.0     	!λͼ
PLNSOL, S,EQV, 0,1.0     	!ʾEQUIVALENT STRESS
ETABLE,P_MX, SMISC, 4    	!аX
ETABLE,P_MY, SMISC, 5    	!аY
ETABLE,P_TX, SMISC, 1    	!аX
ETABLE,P_TY, SMISC, 2    	!аY

PLETAB,P_MX,NOAV         	!ʾаXͼ
PLETAB,P_MY,NOAV         	!ʾаYͼ
PLETAB,P_TX              	!ʾаXͼ
PLETAB,P_TY              	!ʾаYͼ

!ֺ֧
ESEL,S,ENAME,,16		!ѡԪΪPIPE16ĵԪ
PLNSOL, U,SUM, 0,1.0     	!λͼ
ETABLE,Z_FORX, SMISC, 7  	!֧X
ETABLE,Z_FORY, SMISC, 8  	!֧Y
ETABLE,Z_FORZ, SMISC, 9  	!֧Z
ETABLE,Z_MOMX, SMISC, 10 	!֧X
ETABLE,Z_MOMY, SMISC, 11 	!֧Y
ETABLE,Z_MOMZ, SMISC, 12 	!֧Z

PLLS,Z_FORX,Z_FORX,0.5,0        !ʾ֧X
PLLS,Z_FORY,Z_FORY,1,0          !ʾ֧Y
PLLS,Z_FORZ,Z_FORZ,1,0          !ʾ֧Z
PLLS,Z_MOMX,Z_MOMX,1,0          !ʾ֧X
PLLS,Z_MOMY,Z_MOMY,1,0          !ʾ֧Y
PLLS,Z_MOMZ,Z_MOMZ,1,0          !ʾ֧Z

!˲ֺ
ESEL,S,ENAME,,8			!ѡԪΪLINK8ĵԪ
PLNSOL, U,SUM, 0,1.0     	!λͼ
ETABLE,SAXL, LS, 1       	!

PLLS,SAXL,SAXL,0.3,0       	!ʾ

ALLSEL
!-------------------
!--ʱ̺---
!-------------------
FINI
/POST26
NSOL,2,1622,U,X, UX_UP  	!ȡ޶ڵXλ
NSOL,3,1622,U,Y, UY_UP  	!ȡ޶ڵYλ
NSOL,4,1622,U,Z, UZ_UP  	!ȡ޶ڵZλ

NSOL,5,1557,U,X, UX_LOW  	!ȡڵXλ
NSOL,6,1577,U,Y, UY_LOW  	!ȡڵYλ
NSOL,7,1577,U,Z, UZ_LOW  	!ȡڵZλ

DERIV,8,2,1,,VX_UP		!ԶڵXλ΢㣬õXٶ
DERIV,9,3,1,,VY_UP		!ԶڵYλ΢㣬õYٶ
DERIV,10,4,1,,VZ_UP		!ԶڵZλ΢㣬õZٶ

DERIV,11,5,1,,VX_LOW		!ڵXλ΢㣬õXٶ
DERIV,12,6,1,,VY_LOW		!ڵYλ΢㣬õYٶ
DERIV,13,7,1,,VZ_LOW		!ڵZλ΢㣬õZٶ

PLVAR,2,5			!ʾڵXλʱ仯
PLVAR,3,6			!ʾڵYλʱ仯
PLVAR,4,7			!ʾڵZλʱ仯
PLVAR,8,11			!ʾڵXٶʱ仯
PLVAR,9,12			!ʾڵYٶʱ仯
PLVAR,10,13			!ʾڵZٶʱ仯
